Polarons in suspended carbon nanotubes 
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We prove theoretically the possibility of electric-field controlled polaron formation involving fiexu- 
ral (bending) modes in suspended carbon nanotubes. Upon increasing the field, the ground state of 
the system with a single extra electron undergoes a first order phase transition between an extended 
state and a localized polaron state. For a common experimental setup, the threshold electric field 
is only of order ~ 10 -2 V//xm. 

PACS numbers: 73.63.Fg, 71.38.-k, 62.25.-g 



Due to their unique material properies, carbon nan- 
otubes make ideal flexible nano-rods for mechani- 
cal applications^. Coupling their mechanical motion 
to electronic degrees of freedom leads to non-linear 
dynamics 2 . Current technology 3 allows for the fabrica- 
tion of ultra-clean nanotubes in which electrons propa- 
gate ballistically 4 rather than diffusively. In combination 
with a high quality factor 5 , this allows for resonant exci- 
tation and coherent manipulation of discrete degrees of 
freedom. The envisioned devices may find application in 
quantum information processing. 

In current devices, a discrete spectrum is obtained by 
embedding a quantum dot on a suspended nanotube^ 3 -. 
In this paper we prove the possibility of the controlable 
formation of discrete states of a different kind, namely 
polarons. The setup is shown in Figure [T] It consists of 
an ultra-clean carbon nanotube cantilever. We consider 
a single wall semi-conducting nanotube. The setup is 
similar to the nano-relay proposed in Ref. 6 and to the 
experimental setup of Ref. 0, but operated in a different 
regime, namely that of a single electron on the cantilever. 

If the electron enters the suspended part of the tube, it 
experiences a force F = —eE. The electric field E may 
be due to an external source, or to an induced image 
charge in the substrate below the cantilever. The force 
F deforms the tube. As a result, the potential energy 
of the electron is lowered. Thus the tube deformation 
produces a potential well that may trap the electron. 

The trapping of an electron in a lattice deformation in 
a bulk solid is a well-studied topirA The resulting quasi- 
particle is called a polaron. Previous studies of polarons 
in carbon nanotubes^ considered only axial stretching 
and radial breathing modes of the tube, while our study 
concentrates on macroscopic flexural (bending) modes. 

Our main results are contained in Fig. [3] At small 
electric fields, the ground state consists of an undeformed 
tube and an extended electron. As the field is increased 
beyond a critical value, the system undergoes a first order 
phase transition to a localized polaron state. For realistic 
values of a suspended tube length L = 1 jum, and tube 
radius r = 1 nm, the threshold electric field is 0.031 
V//im and the tip deviation is 0.54 nm. This is the field 
that would be produced by an image charge induced in 



FIG. 1: Setup: A single wall carbon nanotube cantilever of 
length L. The supported part of the tube rests on an insulat- 
ing substrate. An electron that enters the suspended part of 
the tube experiences a force F = —eE perpendicular to the 
tube. As a result the tube is deformed so that each point x 
on the tube undergoes a displacement y(x) perpendicular to 
the a;-axis. The electron wave function ip{x) is also indicated. 



a metallic substrate 0.10 /im below the tube. 

We start our analysis by noting that the typical energy 
scale for the polaron state is set by the electron confine- 
ment energy e e = h 2 /2m*L 2 , where m* is the effective 
electron mass. The ratio ?ioJo/e e , where ljq is the fre- 
quency of the lowest flexural mode of the tube, turns out 
to equal 0.09, independent of r or L. This ratio is small 
essentially because electrons weigh much less than car- 
bon atoms. We therefore neglect the zero point motion 
(associated with energy fiu^) of the cantilever and treat 
its displacement as a classical variable. 

The supported part of the tube is tightly clamped to 
the substrate by van der Waals forces and cannot be de- 
formed. We take the tube-element that was at tq = x x 
to be displaced to r — xx + y{x)y under deforma- 
tion. This is valid in the small deflection regime where 
max{|y(x)|} <C L. The system is described by two fields, 
namely, the tube profile y{x) and the wave function ip{ x ) 
of the single electron in the conduction band. The bound- 
ary conditions on the tube profile are y(x < 0) = y'(x < 
0) = and y"(L) = y"'(L) = 0. The boundary condi- 
tions on the wave function are tp(— oo) — tp(L) — 0. The 
wave function can be taken as real, and is normalized. 
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The ground state configuration is obtained by mini- 
mizing the energy functional 
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The term T is the kinetic energy of the electron. The 
effective mass m* is inversely proportional to the ra- 
dius r of the nanotube^. For zig-zag nanotubes m* — 
1.8m e ao/r where m e is the true electron mass and ao is 
the Bohr radius. For tubes with chiralities other than 
zig-zag, the proportionality constant is different, but of 
the same order of magnitude. 

If the electron is at position x in the suspended part of 
the tube, it has undergone a vertical displacement y{x) in 
the direction of the electrostatic force — eE. This means 
that the electron sees a potential well with the same pro- 
file as the tube. The term V in Eq. ([1]) accounts for this. 

The term U is the elastic energy of the deformed 
tube^i. In the small deflection approximation, the en- 
ergy stored in stretching modes is smaller than the energy 
stored in flexural modes by a factor of order [y(L)/L] 2 -C 
1. We therefore only take bending energy into account. 
Y is Young's modulus. It is a material constant, indepen- 
dent of tube dimensions. / ~ irgr 3 is the second moment 
of area of the tube cross-section. Here g is the thick- 
ness of the cylinder wall of the nanotube. Good agree- 
ment with nanotube elasticity experiments is obtained 
by taking g — 6.4 ao (equal to the interlayer distance in 
graphite) and Y = 1.2 X 10 12 Pa-i 2 -. 

It is convenient to introduce dimensionless quantities 
h = ^^H, <f> = VL^J, f = ^y, and z = x/L. The 
dimensionless energy functional h[<f>, f] is explicitly given 
by 



h[<f>J}= f d 
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It depends on a single parameter, the dimensionless cou- 
pling constant a = 2m*(eE) 2 L 5 /h 2 YI. 

Two classes of solutions, or phases, can be distin- 
guished in the system. The first class comprises extended 
electronic states, in which the magnitude of the wave 
function is sizable over the whole length of the tube. (We 
consider a tube with total length 3> L.) For such states, 
the average charge in the suspended part of the tube is 
vanishingly small. As a result, the force exerted on the 
tube by the electric field, and hence the deformation of 
the tube, is zero. The total energy of such a state is equal 
to the kinetic energy of the electron. Therefore the ex- 
tended state spectrum forms a continuum bounded form 
below by zero. The lowest extended state energy is zero, 
corresponding to an electron wave function with an infi- 
nite wavelength. 

The second class of states is of the polaron type. These 
consist of an electron trapped in the potential well asso- 
ciated with the tube deformation that the electron itself 



produces. The electron wave function cf> decays exponen- 
tially into the supported part of the tube, i.e <j> = </>oe K2 
for z < 0, where k is the inverse localization length. Due 
to the negative potential energy of the trapped electron, 
the total energy of the state can become negative. When 
this happens, the ground state of the system is of the po- 
laron variety, since all extended states have positive en- 
ergies. Otherwise the polaron state is meta-stable, since 
there exists an extended state of zero energy. Our task 
is to determine into which of these two classes the groud 
state falls for a given value of a. 

In the limit k<1, where the wave function penetrates 
deep into the susported part of the tube, it is straight- 
forward to estimate the leading order in n contributions 
of the various terms in Eq.[2] (See Appendix |A1 for detail 
of the calculation.) The energy is dominated by a posi- 
tive contribution ~ k of the kinetic energy density in the 
suspended part of the tube. All other contributions are 
~ k 2 or smaller. This implies that the transition to the 
polaron state has to be first order: because lim re _ >0 + 9 K h 
is positive and cannot change sign, the slope of h can 
only vanish at non-zero k. 

Since we are dealing with a first, rather than a second 
order transition, an expansion of the energy in an order 
parameter such as K is of little further use. We therefore 
proceed to a variational calculation. For given <p we find 
the optimal tube profile fo[(f>] by minimizing the energy 
h[(f>, f] with respect to the tube profile /. This is substi- 
tuted back into h to obtain ft. V ar[0] = M^/oML wn ich 
is then varied over a family of trial wave functions (j> vaT . 
We choose a single parameter family 



</>var(z) = N(k) 



e KZ if z < 0; 

[(1 + k)z + 1] (1-z) if < z < 1. 



( 3 ) 

with N(k) 2 — (15 + 16k + 7k 2 + k 3 )/30k ensuring nor- 
malization. The trail wave function has the correct from 
in the supported part of the tube, is smooth at z = 0, 
and satisfies the boundary condition 0(1) = at the sus- 
pended end of the tube. The variational parameter is k. 
Straight-forward but tedious algebra then yields /i var (k) 
as a rational function where both numerator and denom- 
inator are sixth-order polinomials in k. (See Appendix [Bl 
for more detail.) Note as an aside that, to leading order 
in k, we find h e g = 8«/3, consistent with the discussion 
in the previous paragraph. Figure shows /i V ar versus k 
for various values of a. 

The best estimate for the energy for given a is obtained 
by minimizing h vax (Kj with respect to k on the interval 
[0, oo). The following results are found: For a < 350.0 
(dotted curves in Figure ^ , ft, var is a monotonically in- 
creasing function of k so that the minimum is h va-r = 
which occurs at k = 0. The implication is that here 
the ground state is extended. At a > = 350.0, 

a local minimum develops at k > K™[ n = 2.044. For 

a mfn' < a < a c ar = 431.5 (dashed curves in Figure [2]) 
this minimum has a positive energy and therefore cor- 
responds to a meta-stable state. Finally, for a > o>c 



3 




FIG. 2: First order phase transition: The energy /i V ar ver- 
sus the dimensionless inverse localization length k, calculated 
for the trail wave function ([3]), for several values of the cou- 
pling constant a. The dotted curves correspond to a < a^L, 
where the variational calculation predicts no polaron state. 
The upper thick curve corresponds to a = «™ n = 350.0. The 
dashed curves correspond to «™ n < a < a™ where the vari- 
ational calculation predicts a meta-stable polaron state. The 
thick lower curve corresponds to a = «™ r = 431.5. The solid 
curves correspond to a > a™ where the variational curve 
predicts a polaron ground state. Circles indicate the minima 
corresponding to the polaron state. 



(thin solid curves in Figure [5]) , the energy of the polaron 
state becomes negative so that the polaron state is stable. 

Next we numerically minimize h[<f>, f] with respect to <f> 
and /, with appropriate boundary conditions and subject 
to the constraint that <fi is normalized. Details about the 
numerical method can be found in Appendix [Cj 

Thus we find a c = 306.9. (See the top-left panel of Fig- 
ure [3l) This value of a c is lower than the upper bound 
derived by means of the variational calculation above, 
as it should be. It is also of the same order of magni- 
tude as the variational upper bound, indicating that the 
variational calculation is reasonably accurate. We fur- 
ther obtain the value of a m i n , the smallest value of a for 
which polaron states exist as a mm = 254.5 = 0.8293a c . 

At a — a c we obtain a critical tip displacement 
/ c (l) = 0.138. Reinstating units and eliminating the 
electric field in favour of a c we obtain y c (L)/L = 
0.1387i^/ a c /2m*YIL. The critical tip displacement 
scales like r -1 ^ 1 / 2 . For realistic values L = 1 pm and 
r = 1 nm, we find y c {L) = 0.54 nm. 

We also calculate n = U/huJo, where U is the bend- 
ing energy and ujq = 3.52y^Y I / p/ 'L 2 — is the angular fre- 
quency of the lowest harmonic of the suspended tube. 
Here p = 0.671 Mr/a 2 , is the mass per unit length of the 
tube, and M is the mass of a carbon atom. Appendix [Dl 
provides more detail. The quantity n, being the ratio 
between the energy stored in the deformed tube and the 
energy of a single phonon, is an estimate of the num- 
ber of phonons involved in the tube deformation. In the 
lower left panel of Figure |31 n is plotted as a function of 
a. We see that when the transition to the polaron state 
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FIG. 3: Results of the numerical calculation: In the top panels 
and in the bottom left panel, solid curves indicate ground 
state properties. Dashed curves indicate properties of the 
meta-stable polaron state. A thin vertical line indicates the 
value of a/a c below which no polaron solutions exist. The 
critical value of a is a c = 306.9. Top left: The minimal values 
ho of the (dimensionless) energy h[<f>, f] (cf. Eq. ([2}) vs. a/a c . 
Top right: The dimensionless inverse localization length n in 
the suspended part of the tube vs. a/a c . Bottom left: The 
ratio n — U/huuo between the bending energy and the lowest 
phonon energy vs. a/a c . Bottom right: The electric field E 
as a function of a/a c , for a suspended section of length L — 1 
/xm, and three different tube radii r. The solid curve is for 
r = 1 nm, the dashed curve for r = 1.5 nm and the dotted 
curve for r = 2 nm. 



occurs, there are on the order of a hundred phonons in 
the tube. The fact that n is large in the polaron state 
provides additional a posteriori justification for treating 
the tube deformation classically. 
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FIG. 4: Phase diagram: The three different shaded regions 
correspond to different tube radii (assuming zig-zag tubes): 
The dark gray region is for r = 1 nm, the gray region is for 
r = 1.5 nm and the light gray region is for r — 2 nm. In each 
case the shaded region indicates where the ground state is a 
polaron. At large L, the upper planes have been cut away to 
reveal the planes below. At large E the shaded regions fade to 
white to emphasize that the upper boundaries of the shaded 
regions do not represent a phase transition. 

An important question to ask is whether values of a 
that are of the order a c and larger can be reached for 
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realistic values of the length L, radius r, and external 
electric field E. Typical radii are of the order 1 nm. 
Typical lengths are of the order 1 /im. An upper bound 
on the electric field is provided by the breakdown field 
of the insulating elements in the setup. These are typi- 
cally made of SiC>2 for which the breakdown field is ~ 10 
V//im. In the bottom right panel of Figure we plot 
the electric field E versus the corresponding a for L = 1 
/im and three values of r ranging from 1 to 2 nm. We 
see that producing a coupling constant in excess of a c 
requires an electric field of 0.031 V//im for the thinnest 
tubes and 0.12 V//xm for the thickest tubes. These are 
quite reasonable values, well below the breakdown field of 
SiC>2. It is also of the same order as the field produced by 
an image charge in a metallic substrate ~ 0.1 /im below 
the cantelever. 

It is informative to draw a phase diagram, indicat- 
ing the region in parameter space where the ground 
state is of the polaron variety. There are two condi- 
tions that have to be met. Firstly, as we have discussed 
above, the coupling constant must be large enough, i.e. 
a(r, L, E) > a c = 306.9. There is however another con- 
dition: The deflection y(L) of the tube's free tip must be 
much smaller than the length L of the suspended section. 
When \y(l)\ — L, the tube will likely come into contact 
with one of the surrounding elements of the setup, for 
instance the supporting substrate at x = 0. Owing to 
van der Waals forces the tube will adhere to whatever it 
comes in contact with. When this happens, the coupling 
between electron motion and tube profile is destroyed. 
Of course, when this condition is not met, the small de- 
flection approximation, on which our analysis relied, also 
breaks down. Appendix [E] provides more detail. In Fig- 
ure 0] we show three cuts through the phase diagram in 
the E—L plane, for r = 1, 1.5, and 2 nm respectively. In 
each case the polaron phase is indicated by a shaded re- 
gion. We see that the value of the largest allowed electric 
field is always several orders of magnitude larger than the 
smallest allowed electric field. 

In conclusion, we found that the coupling between the 
electron and the tube is controlled by a dimensionless 
coupling constant a — 2m*(eE) 2 L 5 /h 2 YI. At strong 
coupling (large a) the ground state of the system is a 
polaron, i.e. the electron is trapped in the deformation 
of the cantilever that it itself produces. As the coupling 
is decreased beyond the critical value of a c = 306.9, a 
first order phase transition occurs. Below the transition 
point, the ground state consists of an undeformed tube 
and an extended electron wave function. For realistic 
values L = 1 fiia for the length of the suspended tube 
and r = 1 nm for the tube radius, an electric field of 
0.031 V / fxm is required to realize the polaron phase and 
the critical tip displacement is 0.54 nm. The magnitude 
of the threshold electric field is the same as that produced 
by an image charge in a metallic substrate 0.10 /im below 
the cantelever. 

In future work we plan to study the nonlinear dynam- 
ics of a single polaron as well as the interaction between 



polarons in the same or adjacent suspended tubes. The 
eventual aim is to exploit the coupling between mechan- 
ical and electrical degrees of freedom for the coherent 
manipulation of the quantum state of the polaron. 



Appendix A: Expanding the energy in small k 

In the main text it is stated that at small inverse local- 
ization lengths >t, the energy h is dominated by a positive 
contribution of order k. Here we provide a detailed anal- 
ysis. 

Setting 



dz6 2 ) = 



5<p{x) J_ 00 ' ' 8f(x 
and using the boundary conditions 
(f)(z = 1) = lim 4>(z) 

z— > — oo 

/(* = 0) = d z f(z = 0) 



(h-e ( 

J — OC 



o, 
o, 



d 2 f(z = 1) - dlf{z = 1) = 0, 
we obtain two differential equations 



dz(j) 2 ) = 
(Al) 

(A2a) 

(A2b) 
(A2c) 



e(j){z) 
d 4 J(z) 



(A3a) 
(A3b) 



Here £ is a Lagrange multiplier that enforces the nor- 
malization of cf>. The first of these equations is the 
Schrodinger equation for the electron in a potential 
af(z). The second equation describes the balance of the 
electrostatic force that deforms the tube and the elastic 
restoring force. 

The solution to Eq. (|A3bl) that satisfies the boundary 
conditions (|A2bj) and (|A2c|) is 



dz' (z') 2 (3z- z')4>{z'f 



+ / dz' z 2 (3z' - z)4>{z') 



i\2 



(A4) 



We can substitute this solution into h in order to obtain 
an effective energy functional that depends on <f> only, i.e. 
h c g[4>] — h[cf), f,f,]. By exploiting the fact that 



l r i 

2f\2_ / j J f .o4. 



dz'(d 2 u) 2 = dz'udtu = - dz'f^ 2 

o Jo Jo 



we obtain 

hes[(f>] 



dz{dA? + ^ / dz 
I 2 J 



(A5) 



(A6) 



We now expand h c ff in the inverse localization length 
k of the polaron state. Let us firstly look at the kinetic 
energy. We consider separately the contribution of the 
wave function in the supported and suspended parts of 
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the tube. As mentioned before, the wave function is of 
the form </> = (f>oe KZ in the supported part of the nan- 
otube. In the limit k <C 1, the normalization constant (f>o 
is of the order y/H. For the kinetic energy density inte- 
grated over the supported part of the tube, we obtain 




(In the last step we approximated f_ oo dz(j> 2 ~ 1, ne- 
glecting the possibility to find the electron in the region 
< z < 1, which is valid for k < 1.) To estimate the 
kinetic energy stored in the suspended part of the tube, 
we note that here <f> changes from 4>o ~ yfn at z — 
to <f> = at z = 1. This corresponds to a typical slope 
d z (f> ~ — \[k so that 

f dz(d z c/)) 2 ~ K (A8) 
Jo 

Thus, at small k, the kinetic energy is dominated by 
a contribution of order k. Note also that the kinetic 
energy is positive. To estimate the remaining term 
(—a dz fcj,(j) 2 ) in h e g we note from that Eq. (|A4j) that 
f<f,<j) 2 is quartic in <f>. The integral therefore scales like 
<Pq ~ k 2 . It is also negative. 

Appendix B: Variational calculation 
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FIG. 5: Numerically computed polaron states: Tube profiles 
y(x) and wave functions ip{x) that solve Eqs. (|A3a|) and (|A3b|) 
for three different values of the coupling constant a. The 
solid curves are for a = 500 = 1.629a c . Since a > a c for the 
solid curves, they represent the ground state of the system 
at this particular value of a. The dashed curves are for a = 
306.9 = 1.00q c . Since a = a c for these curves, they occur 
at the boundary between stability, where the energy of the 
polaron state is negative, and meta-stability where the energy 
of the polaron state is positive. The dotted curves are for 
a = 254.5 = 0.8293a c = oimm- There are no polaron states 
for a < a m in. 



In the main text we discuss a variational calculation 
in order to obtain an estimate h val (K,) of the energy of 
the polaron state, that depends on a single variational 
parameter k. In Fig. 2 of the main text, h va _ r (n) is plotted 
for several values of the coupling constant a. Here we give 
the explicit formula for /i var («;). It is a rational function 

= EL °?"" Mk " , (bi) 

with coefficients 
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Appendix C: Numerical calculation 

In the main text we present results of a numerical min- 
imization of the energy h. Here we provide some details 
about the numerical method. 



We extremize the energy functional h by solving solv- 
ing Eqs. (|A3a[) and (|A3b[) numerically, subject to the 
boundary conditions (|A2a|) . (|A2bj) and (|A2cp . We use 
an iterative procedure. In each iteration we substitute 
a guess f n {z) for the tube profile into the Schrodinger 
equation (|A3a[) . The associated normalized ground state 
wave function 4> n (z) and electron ground state energy e n 
is computed. The wave function <f> n is then substituted 
into Eq. (|A4[) . This is used as the next guess f n+ i(x) 
for the tube profile, and the process is repeated until the 
change in electron energy \e n +i — £n\ is smaller than the 
required accuracy 5. We choose 5 = 10~ 4 , and an initial 
guess for the tube profile 

A(z) = ~z 2 (3-z). (CI) 

This would have been the exact tube profile (cf. 
Eq. (|A4|) ). had the electron been localized right at the 
tip (z=l) of the tube. In each iteration the total en- 
ergy is also calculated, and we check that it decreases in 
each iteration of the calculation. This gaurantees that 
the obtaind solution is a minimum. 

We find that for a larger than about l.la m i n , conver- 
gence is obtained within less than 20 iterations, while for 
smaller a up to 120 iterations are required. In Figure 
[5] we show (converged) tube profiles and wave functions 
calculated with this procedure for three different values 
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of a. 

The value of a m i n , the smallest value of a for which 
polaron states exist, is obtained by repeating the iter- 
ative numerical calculation for smaller and smaller a, 
until no amount of iteration produces convergence any 
more. We find that for smaller and smaller a down 
to 254.518, the number of iterations required to ob- 
tain convergence slowly increases up to 120. Then for 
a = 254.510, there is a sudden jump, and after 200 it- 
erations, convergence is still not obtained. We conclude 
that a min = 254.51 = 0.8293a c . 



is proportional to the tube radius. (The proportionality 
constant is 1.47 x 10 4 m e / Oq.) Since I is proportional to 
r 3 , wo is proportional to r/L 2 . The bending energy on 
the other hand can be written as 



U 



i,,,>v a J dzU4>2 - 



(D2) 



Since the effective electron mass to* is proportional to 
1/r, m*L 2 u)o is independent of r and L so that n = Uo/uiq 
is a function of a only. 



Appendix D: The parameter dependence of the 
phonon-number n 

In Fig. 3 of the main text the phonon-number n is 
plotted as a function of the coupling constant a. Here 
we show that indeed, n does not depend on L and r sep- 
arately, but only on a, since this is not clear a priori. 
(The above statement holds subject to the approxima- 
tion I = irgr 3 , which is valid when g <C r. A more so- 
phisticated approximation for / yields only a very weak 
r dependence in the regime of realistic radii.) 

Note firstly that ujq is given by 

3 - 52 I™ m-M 

Wo = T^VT' ( } 

where p is the mass per unit length of the tube, which 



Appendix E: Estimating when \y(L) \ ~ L 

In the main text the phase diagram of the system is 
discussed. We limit our discussion of the polaron phase 
to the regime \y(L)\ <~ L. This gives rise to the upper 
boundaries of the shaded regions in Fig. 4 of the main 
text. The estimate for y(L) L was obtained as follows: 
We firstly note that \y(L)\ < |y max | where |y max | is the 
deflection produced when the electron is completely local- 
ized at x = L. From Eq. lCll we have |y max | = eEL 3 /3IY. 
Thus, we demand that 1 > eEL 2 /3IY. 
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